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^ : Abstract 

00 ■ We calculate the density of states of a layered superconductor in which there 

'^ ■ are two layers per unit cell. One of the layers contains a d-wave pairing 

m 

t^ ■ interaction while the other is a normal metal. The goal of this article is 

o. 

to understand how the d-wave behaviour of the system is modified by the 
coupling between the layer-types. This coupling takes the form of coherent, 
single particle tunneling along the c-axis. We find that there are two physi- 
cally different limits of behaviour, which depend on the relative locations of 



Ph. 

^ . the Fermi surfaces of the two layer-types. We also discuss the interference 



between the interlayer coupling and pairing interaction and we find that this 

interference leads to features in the density of states. 
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I. INTRODUCTION 

One of the most interesting developments in the study of high-temperature (HTc) super- 
conductors has been the idea that the gap is not isotropic within the copper-oxide planes. 
There is considerable experimental evidence that this is the case; angle resolved photoe- 
mission experimentsil, for example, have directly observed an anisotropic gap in the aft- 
plane spectrum of Bi2Sr2CaCu208+5. Furthermore, there have been a number of recent 
experiments^ in which the behaviour of the system can be explained by gaps which have 
at least one node on the Fermi surface. There have also been a number of controversial 
tunneling experiments (for example WoUman et afl) in which the sign of the gap has been 
shown to change across the nodes. These experimental results have been taken as strong 
support for a gap with d-wave symmetry within the copper-oxide planes. 

The d-wave symmetry of the gap is an integral part of a number of theories of super- 
conductivity, such as RVBoQ and spin-fiuctuationETEi models. There are also a large number 
of BCS-like calculationsEiEJ (of which we have only mentioned a few) in which the ansatz 
is made that the gap has a d-wave structure. The one feature which is common to all of 
these calculations is that they are 2-dimensional, and it is assumed that the effects of the 
third dimension will be weakly perturbative. This assumption reflects the belief that the 
superconductivity of the HT^ materials occurs within the copper-oxide planes. 

The goal of this work is to consider the effect that the third dimension will have on 
the properties of an otherwise d-wave system. In particular, we will calculate the super- 
conducting density of states (DOS) of a layered S/N system. In this model, there are two 
layers in the unit cell. One of the layers (S) contains a d-wave pairing interaction which 
drives the superconducting transition. The second layer (N) is intrinsically normal in the 
sense that it contains no pairing interaction, although it still becomes superconducting be- 
cause of its proximity to the S layer. Of the HTc materials, this model has the most in 
common with YBa2Cu307, where it is considered that the copper-oxygen planes contain a 
superconducting pairing mechanism and the copper-oxygen chains are intrinsically normal. 



Our particular version of the layered S/N model is based on one studied previouslyOJiia, 
although the focus of these authors was somewhat different than here. In all cases, the 
authors considered the gap of the isolated superconducting plane to be isotropic, and the 
discussion focussed on the anisotropy introduced by the interplane coupling. Abrikosov 
and Klemmll3 were interested in how the interplane coupling affects the density of states, 
the isotropy of the Raman spectrum and the behaviour of Tc and the gap. Buzdin et a^ 
calculated the density of states, while Bulaevskiill3 considered various properties of a layered 
S/N model in the limits of very weak and very strong interlayer coupling. Tachiki et al. have 
calculated a number of observable quantities (tunneling conductance, optical conductivity, 
knight shift and nuclear magnetic relaxation rate) for a more complicated model in which 
there are five layers per unit cell. In this paper, the pairing in the S layer is chosen to have 
a dx^-y'i symmetry and we are interested in how the interplane coupling will change the 
density of states from the 2-dimensional caseuily. 

The layout of the paper is as follows: In section ^ we introduce a model for a layered S/N 
superconductor and derive a formula for the density of states. In section |T^ we calculate 
the DOS for different choices of model parameters. The discussion is broken up into two 
physically different limits. In one case, the Fermi surfaces of the two bands (which result 
from the two layer types) are far apart in the Brillouin zone, while in the second case they 
coincide. The difference between these two limits is that, in the second case, there are 
important (non-perturbative) interference effects between the mean field and the interplane 
coupling. In section |^ we present a summary and discussion of our results. 

II. THE MODEL 

In this section we introduce a model for a layered superconductor, which is based on one 
proposed by Abrikosovtj. In this model we consider two types of metallic layers, stacked 
alternately one above the other in the z direction. These layers form two sublattices which we 
will assume are weakly coupled to each other. There is a BCS-like pairing potential localized 



about the layers of the first sublattice which is responsible for the superconductivity of the 
sample. The Hamiltonian for such a system can be written 

H = E [d\¥,{r){Ho + Vi + \/2)^a(r) 

era' 

where \E'o- is the field operator of spin a. The operator Hq is the free particle Hamiltonian, 
p^ 12771, while Vi and V2 are the potentials associated with each of the sublattices. The 
BCS-like pairing potential is 

/i(r, r') = -1/25: \/(x, y, x\ y')b{z - Z)5{z' - Z), 
z 

where Z are the z-values of the planes in the first sublattice. 

If the electrons are tightly bound to their planes, then the eigenstates of i^o + ^1 + ^2 

are well approximated by eigenstates of the sublattice Hamiltonians Hq + Vi and Hq + V2 

(which are assumed known). Eigenstates of Hq + Vi and Hq + V2 are not orthogonal to each 

other, but may be orthogonalised to form a new set of states, ^jkli"),"^ = 1?2 , which are 

no longer eigenstates of the sublattice Hamiltonians but which satisfy /rf^r 0*j^(r)02k'('^) — 

0. The orthogonalisation procedure will conserve the Bloch-like properties of the original 

eigenstates. If we define 

then the Hamiltonian may be written 

H - N^ = E [alk^aika6(k) + a5kaa2ka6(k) 

kcr 

+ E [^(k)alk,a2k. + r (k)a^kaaika] (2) 

kcr 

- E [^kalktal_ki + Akai_kiaikT_ 



where 



e,(k) = jd?rcP*^{v){H, + V, + V2- /i)0,k(r) 



and 



t(k) = J dh 0tk(r) {Ho + V^ + V2) 02k(r) 



^k = 7T X! ^kk'(ai-k'iaik't). 
" k' 



(3) 



Here ^2 is the volume of the sample and () denotes a thermal average. Equation (|]) is 
the self-consistent equation which follows from making a mean field approximation for the 
pairing interaction in equation (|l|). The form of the potential, Vkk', comes from assuming a 
singlet pairing of electrons of opposite k. We find that 

^kk' = H dxdydx'dy'(t)l^{x, y, Z)(t)\_^{x\ y\ Z) 
z -' 

X V{x, y, x\ y')(l)i-k'{x', y\ Z)0ik'(a;, y, Z), 

which, in the limit in which the electrons are tightly bound to the planes, has no kz or k'^ 
dependence. The interplane hopping term, t, has been chosen to conserve k, which makes 
sense if Vi and V2 have the same periodicity and we have a clean metal. 
In the Nambu formalism, the Hamiltonian can be written 



where 



and 



H - N/i = ^ A^(k)g(k) A(k) + const 

k 



Q 



A(k) 



^IkT 



n-kj 



a2kT 



*2-ki 



a(k) -Ak t(k) 

-Aj; -6(-k) -r(-k) 

r(k) 6(k) 

-t(-k) -6(-k) 



(4) 



(5) 



(6) 



The diagonalisation of the matrix, Q(k), will be simplified if we note that the symmetries 
t(k) = t*(— k) and ^i(k) = ^j(— k) follow from time reversal symmetry. Then the energy 
eigenvalues are Ei = E+, E2 = E^, E3 = — -E_, E4 = — -E+, with 



E. 



2 e.+ci+Ai 



± 



a-a+Ar" 






+ t2[(ei + 6)^ + A: 



(7) 



± 



\ 



e? + el + A^ , ,, 



(t2-66)'-(6Ak)2 



(8) 



(we simplify our notation by understanding t^ to mean |tp and A^ to mean | Akp throughout 
this paper). The first expression for E±, given by equation (|^), is useful for describing the 
weakly coupled (small t) limit, while the second, equation (^, is useful for describing £"_ 
near its minima. 

The Hamiltonian now becomes 



H - N/i = ^ ^ Al{k)A,{k)Ei{k) + const 

k i=l 



(9) 



where Aj(k) is the quasiparticle annihilation operator associated with the i energy band. 
The operator, Aj(k), is defined by 
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A,(k)=$:4(k)A,(k), 

and U{k.) is the 4x4 matrix which diagonalises Q: Uij = Ui{Ej) with 



U.{E, 



{Ej - ^2) A 


-{E,+i,)B 


tA 


tB 



(10) 



(11) 



A = t2_(Ak + i?,+ei)(i?i + 6) 



5 = t2-(A*+i?,-6)(i?,-6 



c = A^[e + {E^ - i^f] + B^[e + {E^ + 6)']. 

We are now in a position to solve for the temperature dependence of the gap, Ak. We 
do this by writing equation (H) in terms of the quasiparticle operators, Aj(k). We get 

^k = ^EE^kk'f/H(k')f/2.(k')/(-i5;.(k')) (12) 

j=l k' 

where f{x) = [1 + exp(/5x)]~^ is the Fermi distribution function. Equation (12) is a self- 



consistent equation which can be solved numerically for the gap as a function of temperature. 
Once the gap at T = is known, then the density of states (DOS) may be calculated. 
The density of states is 

2 



P(w) = -— J d^rIinGf^ir,r,uj) 



where the factor 2 is because we are only counting spin up electrons and N is the number 
of unit cells in the crystal. G^^, is the retarded Green's function: 

/oo 

where 6{t) is the step function and { } are anticommutating brackets. If we write 

2 

i=l 
4 

= 5]{?7i,(k)A,(k,t)0ik(r) 

+f/3.(k)A,(k,t)02k(r)} (13) 

and 



A,(k,t) = e-^^'«*A,(k,0) 



then 



P(^) = ^ E E [Uui^y + U3^i^y] ^i^ - E^m. (14) 

i=l k 

The term proportional to f/ii(k)^ (f/3j(k)^) gives the contribution to the density of states of 
the first (second) sublattice made by the i band. It is therefore possible to write the DOS 
as the sum of the two sublattice densities of states 

p{uj) = pi{uj) + p2{uj). (15) 

We will evaluate both the DOS and the gap equation, (|12D , in the following section. 

III. A DISCUSSION OF NUMERICAL RESULTS 

A. Specification of the Model 

The goal of this section is to discuss the density of states of our model. We will do this 
by examining the features of our quasiparticle energy dispersions. To begin with however, 
we must choose a specific form for our Hamiltonian (equation (|^)). We take 

^i(k) = -2ai[cos{kx) + cos{ky)] - pi, i = 1,2 (16) 

which comes from assuming that the electrons are described by tight binding dispersions 
within the planes. The Brillouin zone is — tt < k^, ky < tt, —-njd < kz < vr/d where d is the 
lattice constant in the z-direction. Since ^i describes the energy of an electron in a Bloch 
state of one of the sublattices, it will have a weak kz dependence. However, a much larger k^ 
dependence will come from t(k), which describes the nearest neighbour hopping of electrons 
between layers of different types. The two chemical potentials in (|1^) allow the two bands 
to be offset from each other. The most significant feature of this choice of ^j is that there are 
saddle points at {k^, ky) = (0, tt) (at which ^j = —pi) and other symmetry related points. A 
great deal has been written about model systems in which the saddle points lie at the Fermi 
surface where they play an important role in the dynamics of the systemtHl. We wish to avoid 
this case since the aim of this paper is to understand the role of the interlayer coupling. For 
all of the results presented in this paper, then, we have taken \pi\, \p2\ > max(|t(k)|, |Ak|). 
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The form of the interlattice couphng follows from the assumption that that the electrons 
are tightly bound to the planes: 

t(k) = tie^^^'^i + tae-*'^^^^ (17) 

Here ti and ^2 are complex constants, and the planes are alternately a distance di and d2 
apart so that d = di + ^2. The condition t(k) = t*(— k) will be satisfied if ti and ^2 have 
the same phase. For our numerical calculations we reduce the number of free parameters by 
setting di = d2 = d/2, ti = ^2 = to/ 2 so that 

t(k) =toCos(M/2). (18) 

One of the things which will become apparent is that, although the fine structure of certain 
features of the DOS depend on the specific form we choose for t(k), the existence and 
important properties of these features are independent of this choice. 

The pairing potential is taken to be separable so that \4;k' = Vv^Vk'^ with //k = 1 for 
s-wave superconductors and ?7k = cos(A;^) — cos{ky) for d-wave superconductors. In this case 
we can define Aq by 

^0 = ■^I]^k(ai-kiaikT), (19) 



and equation (|T2D will simplify to 



^0 = ^EE^kf/i.(k)t/2.(k)/(-E,(k)), (20) 

i=l k 

with Ak = Agr/k. We solve equation ( pOD for Tc and graph the results in Fig. |[ We can see 
that the critical temperature decreases rapidly with increasing to which suggests that a weak 
coupling of the planes is necessary for any reasonable description of a HT^ superconductor. It 
has been suggested"'^, however, that the inclusion of phonon mediated interlayer coupling 
may actually lead to an increase in Tc with coupling strength. We also calculate Ao(T = 0) 
from equation ( pil| ) and plot Aq(T = 0)/Tc in Fig. |I|. We can see that Aq(T = 0)/Tc 
increases rapidly for to larger than Aai. For YBCO the value of 2Aniax/^c depends on the 



experiment from which it is determined, but it is generalljH found to be less than 8. It is 
important to remember that the experimentally measured value of the gap is typically the 
lowest energy peak in the quasiparticle excitation spectrum, which is not necessarily related 
in any simple fashion to the order parameter. However, as a rough estimate, we can limit 
max(2Ak) = 4Ao < 8 so that, from Fig. |I|, t^/ai < .5. 

Once Ao(T = 0) is determined, we can proceed to calculate the DOS for our Hamiltonian. 
Our discussion of the DOS will be aided by two observations. The first is that equation (p!6D 
implies that 

6 = ^ei + ^^i^^^i^ (21) 

CTi (Ji 

and the second is that the Brillouin zone boundary, kx,ky = ±7r maps into the diamond 
shaped boundary 



2Af 



1 



(22) 



with — 4crj — Hi < ^i < Aai — /x, and i = 1 or 2. Equations ( pi| ) and ( P^ allow us to discuss 



the system using ^i (or ^2 for that matter) and A^, rather than k^ and ky, as independent 
variables. 

We will consider two limiting cases of Eqn. (pT]) . In the first case |/iicr2 — At20"i|/o'i ^ Aq 
so that ^f + ^1 ^ A^ everywhere in the Brillouin zone. In the second case, we take the 
Fermi surfaces of the sublattices to coincide, so that /iiO"2 ■" Ai2<7i = 0. The fundamental 
difference between the two cases is that it is possible to distinguish band structure effects 
from superconductivity effects in the first case, while in the second case the interplane 
coupling plays an important role at the Fermi surface. 

B. Case I: Distinct Fermi Surfaces 

We begin by discussing the quasiparticle energy dispersions. First of all, for most regions 
of the Brillouin zone ^l + A^ — ^| ^ t^, A^ and we see from Eqn. (|^ that E± approaches 
the limiting form 
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E+ ~ max( J^i + A2, |6 



E_^mm{Jei+Al\^2\). (23) 



In Fig. 1^ we plot E± along the line k^ = ky, kz = and compare it with |^2| and J^f + A^. 
We have chosen an s-wave gap for illustrative purposes in the figure, since a d-wave gap 
would vanishe along this line. In regions of the Brillouin zone where (^) is approximately 



true, the effect of the matrix elements Uij will be to associate JC,j + A^ with the first 
sublattice and |^2| with the second sublattice, so that the integrands in equation (|I^ ) will 
become 



1 + 



^M^ 



6{uj - J a + A2) 



1- 



6 



ViFA 



5{u + v'e? + A^) (24) 

for pi and 

5{u - 6) (25) 

for p2- For most energies, then, the DOS will simply look like the DOS of the decoupledSii3 
(small t) limit, which we show in Fig. |^ 

It is the regions where the interlayer coupling has an important effect which we will 
discuss in the most detail. In Figs. ^ |^ and ^ we show the densities of states for nonzero 
interlayer coupling. These have a number of features that are absent in the uncoupled case, 
and we will discuss these features in turn. 

The first thing which is clear is that, as to is increased and the system goes from be- 
ing 2-dimensional to 3-dimensional, the logarithmic divergences associated with the van 
Hove singularities (both the superconducting one and those intrinsic to ^i and ^2) become 
broadened peaks^. This process is particularly clear in Figs. |^ and ^ where the intrinsic 
singularities at u; = .8 in pi and u = —.4 in p2 are broadened considerably by to- 

The interlayer coupling is also important in regions of the Brillouin zone where t^ is of 
the order of C,f + A^ — ^|. The effect of the coupling is to create avoided band crossings 
wherever 
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^l + Al-e2 = 0- (26) 

In Fig. 1^, for example, we can see that there are two avoided crossings, one at k^ = ky = 1.57, 
and one at k^ = ky = 1.32. 

At one of the avoided crossings {kx = ky = 1.57 in Fig. H), ,^i and ,^2 have the same sign 
and (expanding E± in powers of Aq) 

i?±~|e'|±|t(k)|+0(|l), (27) 

where ^1 = ^' is the solution to Eqns. (|26|) and (0) for which ^1 and ^2 have the same sign. 
The energy difference of the bands at the avoided crossing is ~ 2|t(k)|. We can see from 
Fig. ^that this avoided crossing does not introduce a van Hove singularity and, as a result, 
we do not find any feature in the DOS of Figs. ^ and |^ at this energy. On the other hand, in 
Fig. 1^, we find that, because we have let /ii = /i2, the avoided crossing occurs at the intrinsic 
van Hove singularities of ^1 and ^2- The result of this is that the logarithmic peaks evident 
in Fig. ^ at LJ ~ .8 are reduced far more than in Fig. ^ where the avoided crossing is far 
from the singularities. This is an important result because it shows that models in which 
the van Hove singularity plays an important role (ie. it is near the Fermi surface) cannot 
also have a band crossing near the Fermi surface. 

For the second avoided band crossing (at kx = ky = 1.32 in Fig. (H)), ^1 and ^2 have 
opposite sign. Again we can expand E± in powers of Aq, and we find that the avoided 
crossing occurs at 



2 



^*-Ve" + 'M^ ±7^1)5 + «'!?' (2^) 



(^1 = ^" is the solution to ( p6D and ( pTD for which ^1 and ,^2 have opposite sign). We can see 
from Fig. |^ that there is actually a band gap of width 

E,.E_^^m^,0§) (29) 



at this point. What is interesting about this gap is that it only exists in the superconducting 
state (since it is proportional to A^) and that it forms away from the Fermi surface. This is 
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unusual since the effects of superconductivity normally manifest themselves near the Fermi 
surface. In Fig. ^ these structures are clearly visible at cu ~ S," = 0.55 in the first sublattice 
and uj ~ —C," in the second sublattice. In Fig. |^, the larger value of to shifts the location of 
the minima slightly to u; = ±0.65. In Fig. ^, the structures are clearly visible at tu = ±0.4. 
The widths of the gap-like structures are given approximately by 



E, 



E_ 



4tnA, 



0^0 



max 



e" + ti 



1 



1/^1 ±ri 

4cri 



(30) 



since the maximum value of Ak allowed for ^i = ^" is given by Eqn. (p2D. The source of this 
gap is not immediately obvious, but the proportionality of Eqn. (^) to toAo suggests that 
it results from an interference of interlayer coupling and the mean field. In section [IIID 



we 



will show that this is in fact the case, and we will discuss why these gaps appear on opposite 
sides of the Fermi surface in each of the sublattices. 

The final place where the interlayer coupling is important is near the Fermi surface, 
where it influences the superconducting properties of the system. In Fig. ^ we see that E_ 
has two local minima which are related to the superconducting gaps of the two sublattices. 
We expand E_ (Eqn. (|)) in (t^ - ^^^2) and A^ to get 

it' - 66)^ ± (Ak6)^ 



E^ ~ 



a + ii + 2t^ 



(31) 



Equation (|31| ) will vanish whenever Ak = and t(k)^ = ^1^2- This happens for two values 
of ^j, which we denote by ^^ with j = 1,2. From Eqn. (^) it is straightforward to show 
that 

0) _ /i2tTl - Hia2 



Ci 



2ao 



+i-r 



Aj) 

S2 



+ i-l' 



J V^(/i2fXl - 


-/il(T2)2±4ai(T2t2(k) 


-/UlO-2 
2(Ti 


2^2 


. 7(/i2CXl - 


-/iia2)2±4ai(T2t2(k) 



(32a) 



2ai 



(32b) 



We then estimate the sizes of the superconducting gaps by substituting the equations 
^i = ^p'^ and t = to into Eqn. (|3l|): 
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E_ 



2An 



f(i) 
Aj) _i_ AJ) 

SI "I" S2 



4cri 



(i)i 



(33) 



In the limit that to becomes small, the gap sizes become 



E_ 



2An 



1 



\M 



4cri 
Aot(k)2 



2(/i20"l - /iicr2) 



sgn(/ii) - 



and 



E_ ~ 2Ao 



0"2to 



2 r 



(7l(4eri - |/ii|) 
fl2Cri - yUiO-2 



1/^2 I 

4ao 



(34) 



(35) 



_/i20-i -/ilCr2_ 

Equation (|3^) shows how a weak interlattice coupling affects the gap in the intrinsically 
superconducting sublattice. One interesting point is that whether to increases or decreases 
the gap size depends on the details of the band structure. When to = we have an approx- 
imate expression for the location of a logarithmic singularity in a two-dimensional d-wave 
superconductor. The exact expressions is Ao(4cri — IfiiD/JAaf + Aq. Equation (|35| ) shows 
the size of the induced gap in the intrinsically normal sublattice. 

The gaps described by equations ( [5^ ) and ( P^D will appear in both pi and p2 (this has also 
been pointed outEJ for a similar model). The magnitudes of the contributions are determined 



by the matrix elements, f/jj, in Eqn. ([1^) . In the appendix, we derive the low energy DOS 



and can see explicitely how the matrix elements behave. In particular, Eqns. ( [A9D describe 
the strength with which each of the minima of E_ (labelled by i = 1, 2) is reflected in pi 
(from Uf) and p2 (from f/|). In the weakly coupled (small to) limit, for example, equations 
(^) tell us that f/f ~ 1, f/| ~ (o"it)^/(/i20"i — /ii(j2)^ for the intrinsic gap given by (p4D, 
while Uf ~ (o"2^)^/(Ai20"i — /^iO"2)^, ?7| ~ 1, for the induced gap given by (|35D. In other 
words, features associated with one sublattice are reflected in the other sublattice with a 
strength proportional to tg. We can see quite clearly in Fig. |^, for example, that both gaps 
are reflected in both densities of states. 

The main results of the appendix are expressions for the low energy DOS in which a) 
two of the three k-space integrations have been performed approximately (Eqn. |A10| ) and 
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b) the final integration is performed in a crude approximation in tlie small to limit to find 
Pi (Eqns. ([A16| )) and p2 (Eqns. ( [Al?] )). The important features of these expressions can be 
summarised simply. First of all, the quantity t*{uj) (Eqn. ( [A15| )) gives the energy scale of the 
induced gap. In fact, the equation t{kz) = t*{uj) is just a rearrangement of the expression 
(P5|) for the induced gap. Eqn. ( |A16aD (with X± given by Eqns. ( |AllaD and ( |Allb| )) is an 



approximate expression for the d-wave gap in the S sublattice. Eqns. (|A16b|) and (|A16c|) are 
the corrections (of order tg) to pi inside and outside the induced gap respectively. Similarly, 
Eqns. ( |Al7a|) and ( |A17b|) give p2 inside and outside the induced gap respectively. One of 



the most important features of the induced gap is that, since both t{kz) and Ak have nodes, 
there are a large number of low energy excitations. The density of states still vanishes inside 
the induced gap but instead of vanishing linearly, it vanishes with a divergent slope. We can 
see this in Fig. |^ where we compare the low u DOS, calculated numerically from Eqn. (^4|), 
with the approximate expressions derived in the appendix. 

C. Case II: Coincident Fermi Surfaces 

We will now consider the case where the Fermi surfaces of the two sublattices coincide, 
which requires that piO"2 = P20"i so that Eqn. ( pi]) becomes 

6 = ^. (36) 

Away from the Fermi surface, the structure of the DOS is that of the uncoupled (to = 0) 
limit for exactly the same reasons as in the previous case. Near the Fermi surface, however, 
the DOS behaves quite differently and is somewhat more complicated to describe, though a 
partially quantitative understanding is still possible. In Fig. || we can see that E^ still has 
a double minimum structure throughout much of the Brillouin zone but that it is no longer 
possible to associate either of the minima with a particular sublattice. In Figs. |^ and [1^, 
we show the DOS for two different strengths of interlayer coupling. In Fig. ^, the density of 
states exhibits a nested gap structure at the Fermi surface. For larger values of to (Fig. jTU]), 
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the outer structure is no longer identifiable as a gap and the DOS is similar in appearance 
to that of a 2-diniensional d-wave material. However, as we shall see, the dependence of the 
gap width on the order parameter, Aq, is vastly different than in the 2-dimensional case. 
The complicated structure of Figs. ^ and |1^ is due to the mixing of band structure and 
superconductivity effects at the Fermi surface. Because of this, it is not possible to replicate 
the analysis of the last section where we were able to identify gap-like structures in the DOS 
with features of the energy spectrum. Instead, it is necessary to resort to the brute-force 
method of identifying saddle point singularities in E±, and determining which features in 
the figures they are responsible for. The method used below is approximate and works well 
provided that Aq/cti ^ 1. 

Since finding the zeros of Vk-^i is difficult, our approach is to set dE±/d^i = to find 
the surfaces of extrema of E± with respect to ^i, and then to assume that the van Hove 
singularities are at the extremal values of |t(k) | and | Ak| on this surface. The extremal values 
of |t| are simply zero and Iq, which is actually what we would find by setting Vk-E± = 
anyway. On the other hand, the extremal values of | Ak| are either zero, or given by Eqn. (P^. 
Since this last constraint depends on ^i, it is wrong to treat Ak as independent when taking 
dE±/d^i. Providing that Aq is small, however, the approximation is good and introduces 
an error of 0(Aq/(t^) to the energies of the van Hove singularities. 

We begin by substituting equation (|36D into the expression for E± (equation (^) and 



setting dE±/d^i = 0. We find that there are two solutions: 

^1 = (37) 

(for both E+ and -E"-), and 



_ aijaf + alMk)\^t{ky{a, + a^f + 2a2A^(ai - a^) - ola2[AUai - a^) + 2t(k)^(ai + a^)] 

(72 (CTI +0-2) (di -0-2)2 

(38) 

(for E- only). This second equation gives the location of the two minima of E^ and only 
has positive solutions for ^\ if 
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tikf > 



a^Ai 



((71 + a2)iaf + ct|) - 2af + {aj + (yl)Jal + 2aia2 + 2af 



(39) 



{(^l + 0-2 J 

The important difference, then, between Figs, ^(a) and ^(b) is the size of |t(k)|. This 
crossover in behaviour was noted previouslytj, for the slightly simpler model in which ^i = 1^2 

and Ak = Aq. Now at ^1 = 0, 

Ak 



E+ = \ [-^ 



Ak^^ 



t^± 



(40) 



while when ^1 is given by (^ 



E^ 



2|t(k)|(TiCT2\/t(k)2((Ti + cr2)2 + 2a2A2((Ti - (Ta) - CT2[A^ff2(ffi - (X2) + 2t(k)2ai(ai + ^2)] 



(0-1 +0-2)(cri -0-2)^ 



(41) 



We first describe the saddle points at which ^1 = 0. The extremal values of |t| are |t| = 
and |t| = to) and the extremal values of Ak are zero and (from (|22D) 



A' = 2Ao 



|/ii| 
4cri 



(42) 



Of the four points in the Brillouin zone just described, only two are actually saddle points 
(as opposed to band minima). There will be a saddle point in Ej^ at Ak = A', t(k) = at 
which 



E , = A'. 



(43) 



There will also be saddle points in Ej^ and E^ at Ak = A', t(k) = to at which 

A' 



E. 
E, 



y 



/2 
'^0 



^) +tl + ^. 



2 

A^ 

y 



(44) 
(45) 



The double minimum, with energy given by (^T]) will also have two saddle points, though 
they are somewhat more complicated to describe. First of all, there will only be saddle 
points if the inequality (^) is satisfied for t = to aiid A = A". The quantity. A" is the 
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solution to the pair of equations, (|3^) and (p2[) with t(k) = to- These equations can be 
solved iteratively in a straightforward fashion, and will yield two closely spaced solutions 
corresponding to the two slightly different energies at which the minima of E^ intersect the 
Brillouin zone boundary. The energy of the saddle point is then given by Eqn. (p|). 

In Fig. y. A' = 0.14, and the singularities at ±A' can be associated with the outer gap 
in the nested gap structure. If we recall Eqn. (l3^ ) and the discussion which follows, we 
will recognize that A' is approximately the energy of the gap in a 2-dimensional d-wave 
superconductor. In Fig. |10|, however, the nested gap structure is no longer clear and, and 
although there are features at ±A' = ±0.23, they do not appear as a gap-like structure. 
In Fig. |1^, the gap width is actually given by Eqn. (^T]) which gives the closely spaced sets 
of singularities u = ±0.072 and u = ±0.096. Eqn. (|4l|) makes it clear that the structure 
in the DOS which we naively associate with the superconducting gap is, in fact, due to a 
complicated mixing of the gap and the interlayer coupling. It is an important point that 
a tunneling experiment might not be able to distinguish between Fig. |10| and the density 
of states of a 2-dimensional d-wave superconductor (pi in Fig. |^). In this case, then, a 
tunneling experiment would reveal very little about either the interlayer coupling or the 
nature of the pairing interaction. 

In Fig. ^ the singularities due to Eqn. (|4l|) are at u = ±0.0431 and u = ±0.0426. They 
are responsible for the inner gap in the nested gap structure. There is also a van Hove 
singularity nearby, at u = 0.05, which comes from Eqn. (^^. It is difficult to see because of 
its close proximity to the other singularities, but in Fig. ^, this feature is clear and is at cj = 
±0.30. Once again, this feature is the result of the interference of interplane coupling with 
the mean field and it cannot be considered a perturbation of either a superconducting gap or 
a band gap. In Fig. ^ Eqn. ( ^Sj) describes the smearing of the logarithmic singularity due to 
the intrinsic superconductivity of the first sublattice. It gives van Hove singularities at u = 
±0.195. This should be compared with the second term of Eqn. (p4D, whose kz dependence 
also smears the logarithmic singularity associated with the intrinsic gap. The difference is 
that, in (|34D , t appears as a perturbation, while in (|45|), the smearing is inherently non- 



perturbative. In Fig. [1^, the van Hove singularities are at u; = ±.534. The large broadening 
of the logarithmic singularity has caused it to overlap adjacent features so that it is not 
obvious that it is related to the singularity at A'. The remaining features in the DOS of 



Figs. and 10 are due to the inherent structure of ^i and ^2- 



D. Discussion of the Superconducting Band Gap 



In section [1IIB| we discussed the sublattice densities of states in the limit |/iicr2 — /i2cri| 3> 
AqCTi. We found that, in addition to the expected band and superconducting gaps, there 
were gap-like structures which appeared in the superconducting state away from the Fermi 
surface. This is surprising since the effects of superconductivity normally manifest them- 
selves at the Fermi surface. A further property of these structures is that, if in one sublattice 
a gap-like structure should appear at some energy, u, then there will be a similar structure 
in the other sublattice at energy —u. The purpose of this section is to understand the source 
of these unusual gap-like structures. 

The first clue to the nature of the structures is in the expression (^), which shows that 
the gap width is proportional to Aoto- This suggests that there is an interference between 
the interlayer coupling and the mean field. In order to examine this interference further, we 
will write out an equation for the Green's function in which the interlayer coupling and the 
mean field are treated as perturbations. The Hamiltonian, Eqn. (|), can be written as 



H - N/i = ^ A\k){Q° + Q^}A(k) + const 



(46) 



where 



and 



Q' 



ei(k) 

-ei(-k) 

6(k) 

-6(-k) 



(47) 



19 



Q' 



(48) 



-Ak t(k) 
-A^ -r(-k) 

r(k) 

-t(-k) 

Equation (^) is treated like a perturbation. The temperature Green's function is defined 

as Gij(k; r) = — (TAj(k, —ir), A](k, 0)), where T is the time-ordered product and the time 
dependence of Aj(k, t) is determined by the full Hamiltonian, Q. It is straightforward to 
show that 



G,Mn) = GlitCn) + E GlitCn)QLGmMn)- 



(49) 



l,m=l 



where G{i(n) and G^{i(n) are the Fourier components of G{t) and G^{t) and (n = 'n'{'2n + 
l)//3 are the Matsubara frequencies. The uncoupled Green's function, G^-, is the full Green's 



01-1 



function with to = Aq = 0. It is easy to show that G^j{i(n) = Sij[iCn — Q 

With a little work we are able to write out equation ( ^91) for Gn explicitely: 



Gil — Gil + GiiA]s_G22^l_Gii + GiitG^^t Gu 
+GiiA]f_G22t*G4i 



(50a) 



j-r^o +*/ 



G41 — G44tG22Aj,Gii + G44tG22^ G41. 



(50b) 



The Green's functions, Gn and G33, describes the propagation of spin-up electrons in the 
first and second sublattices respectively. On the other hand, G22 and G44 describe the 
propagation of spin-down holes in the first and second sublattices. The unusual term, G41, 
is essentially defined by Eqn. ( ^Ob|) and it describes an electron which both hops between 
the planes and interacts with the mean field. 

In Figs. |ll|(a) and |ll|(b) we write Eqns. ( 50a ) and ( |50b|) out in diagrammatic form. In 
part (a) we see that there are three processes which modify the propagation of an electron. 
The first of these processes, which simply describes the hopping of electrons between adjacent 
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layers, would appear in any band structure calculation involving two atoms per unit cell. 
The effect of this term is to cause the two bands, .^i and C,2, to repel each other wherever 
they cross so that we have avoided band crossings whenever .^i = ^2- The second process 
in Fig. 0(a) describes the intrinsic superconductivity of the first sublattice. It is useful, 
however to interpret this in the same way we interpreted the first process: as a coupling of 
two bands. The important point here, however, is that the two bands which are coupled 
by the mean field are not the two sublattice bands, but are, instead, the spin-up and spin- 
down components of the band in the first sublattice. Moreover, the coupling is actually 
between spin-up electron states and spin-down hole states. The terms in the Hamiltonian 
responsible for the superconductivity can be interpreted as describing transitions between 
the electron and hole bands in exactly the same way that the terms responsible for interlayer 
coupling describe transitions between electron bands. This is clear in Fig. |ll|(a) where an 
interaction with the mean field changes an electron into a hole. The second term in Fig. 
|Tl](a), then, will cause an avoided crossing in the quasiparticle energies wherever the original 
electron and hole bands cross. Since the dispersion of the original hole band is just — ^i, the 
superconducting band gap will open up at ^i = — ^i = 0. We can see that the reason that 
the effects of superconductivity normally manifest themselves at the Fermi surface is that 
the interaction between the hole and electron bands is largest there. 

The final term in Fig. 0(a) is interesting because it describes the interference of the 
mean field and interlayer coupling. In (c) we have iterated this once, using Fig. 0(b), to 
show the lowest order contribution to the mixing. The diagram in (c) describes an electron of 
energy ^i which interacts with the mean field to become a hole of energy — ^i and then hops 
to the second sublattice before reversing the process. The effect of this diagram is simply to 
couple electron states in the first sublattice with hole states in the second sublattice. Exactly 
as before then, we should expect avoided crossings in our quasiparticle energy dispersions 
whenever ^i = —^2- As we discussed in section p.IIB|, this avoided crossing is responsible for 



the gap-like structures which appear away from the Fermi surface in Figs. ^, ^ and ^. 

In section |IlI Q we discussed the case of coincident sublattice Fermi surfaces. In this case. 
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Fig. |TT|(a) still describes the system, but all three processes will contribute to the avoided 
band crossings at the same time. In particular, the final process in (a) will be responsible 
for the complicated nature of the gap at the Fermi surface. 

IV. CONCLUSION 

In this paper we have presented results on the density of states of a layered S/N sys- 
tem. The inherently superconducting sublattice was presumed to have a gap with d-wave 
structure. We investigated how the (coherent) coupling between the two types of layers 
changed the DOS from the usual 2-dimensional case. Our chosen Hamiltonian, Eqn. @, 
led to a relatively complicated quasiparticle energy dispersion (Eqn. (|^ or (H)), which we 
examined in an attempt to understand our numerical results for the DOS. We found that 
the behaviour of our model could be summarised by two limiting cases, distinguished by the 
relative positions of the sublattice Fermi surfaces. 

In the first case, the Fermi surfaces of the S and N sublattices were far enough apart in the 
Brillouin zone that the interlayer coupling did not cause them to interfere with each other. 
In this limit it was clear that the effect of the interlayer coupling was to shift the usual 2- 
dimensional d-wave gap in the S sublattice perturbatively (Eqn. (0)) and to induce a small 
gap at the Fermi surface of the N sublattice (Eqn. (|35|) ). The gaps in each of the sublattices 
were also weakly refiected in the DOS of the other sublattice, with a magnitude proportional 
to tg. In the appendix, we derived an approximate expression for the density of states within 
the induced gap (Eqn. |A17a|) , and we found that the behaviour is much different from the 



linear behaviour found in the unperturbed d-wave case. The most interesting features, 
though, of the DOS in this case were the gap-like structures which appeared away from the 
Fermi surface in the superconducting state (this was clearly visible in Fig. ^ and we gave 
an expression for the gap widths in Eqn. (|29D). In section [1I1D| , we discussed this matter 



in some detail. Essentially, the point was that there is an interference of the interplane 
coupling and the mean field. Since the mean field can be viewed as an interaction between 
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holes and electrons within the superconducting sublattice, and the interplane coupling is 
an interaction between electrons in the two sublattices, the interference of the terms results 
in an interaction of holes and electrons in different sublattices. It is this interaction which 
resulted in the gap-like structures. 

In the second case we considered, the sublattice bands (^i and ^2) crossed at their Fermi 
surfaces. In this case, there was a nontrivial mixing of mean field and superconductivity 
effects at the Fermi surface. In Fig. ^, we showed that, for weak interlayer coupling, the 
system has a nested gap structure, where the outer gap is associated with the usual d-wave 
gap, and the innner gap is a complicated function of the coupling and order parameter (Eqn. 
(PI)). As the coupling was increased, however (Fig. |10D, the outer gap became unrecognizable 



and the total DOS became very similar to that for a single layer system in which there is 
only one gap. 

Perhaps the most important generalization which follows from this paper, then, is that 
the interference of the interlayer coupling and the superconducting mean field can produce 
surprising behaviour. One implication of this is that for a system like YBa2Cu307, where 
gap-like features have been observed in tunneling and Raman scattering experiments, it 
could be very wrong to simply associate these features directly with the order parameter. 
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APPENDIX: LOW ENERGY DOS 

In this section we will derive the low energy approximation for the density of states in 
the case in which \ai^2 — o"2Aii|/(o"i + ^"2) ^ ^k, ^- One immediate consequence of this 
condition is that we can ignore terms in the DOS (equation (|14|)) containing 6{w±E^) since 
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Ei>ma^{ei,e2)> 



(Ti/i2 - CT2/il 



(Al) 



<Tl + <T2 

This leaves us to consider the contributions to the DOS made by E^. Eqn. ( pl| ) is our first 



approximation to E^ at low energies and we recall that it is zero whenever ^i = ^} , ^2 = ^2 > 



and Ak = with i = 1,2 and Q-^ given by Eqn. 



We will proceed by considering the 



contributions from the neighbourhoods of the two zeros of E_ separately, and we break the 
integral in Eqn. ([l^) into a sum of two integrals which are centred about each of these zeros. 
We define the coordinates 



r/« = 6 - ^P 



and expand E_ to lowest order in r]^^^ and A^: 



El ~ a«2 (.)2 ^ ^(i)2^2 



(A2) 



(A3) 



with 



(.) _ O-id') + ^2^P 



6« 



Ji)2 
S2 



t2+#^' 



The Jacobian for the transformation (k^, ky) -^ {f]^'^\ A) (where A = Ak) is 

1 



J 



4(71 Ao sin(A;^) sm.{ky) 
Ao 



aiJ[«^-A2][«^-A2] 



(A4) 



with 



oi± = 2Ao 



1± 



4(71 



Wi 



This transformation maps the square < kk,ky < tt into a diamond shaped region whose 
boundaries are given by the lines A = ±q;_, — 4cri — /ii — ^| < t]*^*) < Aai — yUi — ^| . Then 



QiK'T / CI thy 





t(») 



f(») 



4<7i-Mi-Cl' 
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(ir/ 



(i) 



dAJ(r/W,A). 



(A5) 



Since we are interested in small tu, we make the further approximation that 



a± ~ 2Ac 



1± 



l/ii + e 
4cri 



Wi 



(A6) 



We can now interchange the order of integration in (|A5|) so that for small positive cu 
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w/d 



pM=N ^ Y. Ul,-i)2^{w - E_] 



2AnrfA ri'^ 



aiTT-' 



E 

2=1 



dk^ 



dA f ""'''' ^.W^li^^ 



dr]^ 



U'.,6{w - E_) 



-4-i-Mi J[al-/\'^][al-/S? 



(A7) 



Finally, in the limit of small uo, 

Ul2 



2 1 ^2 



(02 



f/^ 



1 t" 

32^2t2 + ^W2 



a»ry(i) 



1 + 






UJ 



(A8a) 
(A8b) 



The term a^''^ r]^'^'^ / cu is not small but it is odd in r/^*^ and will vanish in the integration. 
We perform the integral over i]^^^ in (^) to get 

dA 



Pj[UJ) 



?^±rik, "' '^- 



criTT-^ 



i=l 



iW6«| 



[a% - A^][al - A^E^y - A^] 



where Amax = min(a_,ci;/6*^*^) and 



U', 



u' 



M)2 
S2 






(A9a) 
(A9b) 



If Amax = oj/U-"^' then the energy surface E^ = cj is completely contained within the Brillouin 
zone. The integral over A can also be done analytically to give 



Pj\ 



UJ] 



*j^±r,Lii^K 



cTivr'" 



j=i 



i«ka: 



A, 



(AlO) 



The function, K{x), is the complete elliptic integralEll. If cu < a^b^^^ then 

A± = a^[(a_6«)2 + u^] - 2{uja_)^ 



(Alia) 
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while ii LJ > aj-fe*-*) then 



±2(a+a_)y'[cu2 - (a+6»)2][tu2 - (a_6W)2]. 



(Allb) 



There is nothing special about cu = a+fe*^*-* and the integrand in ([A10|) is featureless at 
these points. On the other hand, A_ vanishes at cu = a^b^^^ and K[{\^ — A_)/A+] diverges 
logarithmically. The integration over kz reduces the divergence to a peaked structure. The 
magnitude of these peaks in each of the densities of states is is dictated by Uf. 



The complexity of ( |A10| ) prevents us from performing the final integral analytically 
without further assumptions. So far we have used the fact that Aq and uj are small and we 
now add to this the assumption that to is small. For one value of i, 



;-l)*sgn(/i2cri - /Uicr2) = 1, 



(A12) 



and ^i^ ~ (/i20"i - /iiO"2)/(T2, ^2 ~ o"2i^/(yU20"i - 1^1(^2), a^*-* ~ o"2/o"i, b'-'^ = Uf ~ 
((T2t)^/(/i20"i — /^i'72)^, U2 ~ 1, and a± ~ 2Ao[l ± |/i2|/4a2]. For the remaining value 
of i, 



(-l)*sgn(/i2ai - /iia2) = -1, 



(A13) 



and C,[ 



(i) 



(Jit^/{l22(Ti - /iia2), G {1^2(^1 - /ii^2)/f^i, a^'^ ~ b^^ = Ul ~ 1, f/| 



(crit)2/(/i2cri — fiio-2Y, and a± ~ 2Ao[l ± |/ii|/4cri]. Substituting this into (|A10| ) gives 

"A^-A. 



4Aocj 1 ^^ 



A, 



+ 



AdAouj (T2 

7r3 [/i2cri - /xicr2]2 Jo 






A+-A_ 



(AM) 



The first term is for i given by ( [A13| ) and has \± given by ([Allap since 6^*^ ~ 1. The second 
term is for i given by ( [A12| ) and, since b^'^' is a function of kz^ A± will be given by ( |Allaj) 
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when t{kz) > t* and by ( |Allb|) when t{kz) < t*. The value of t* conies directly from the 
condition a; = a^b^^^: 



t*{u;) 



2uj 



lyU-20'l -/iicr2|. 



(A15) 



/ AoO-2[4o-2 + \lJ,2\] 

In order to solve the final integral in ( |A14D , we make the crude approximation that X± can 
be replaced by the limiting forms X±(t ^ t*) = cj^[a+ ±a_]^ when t < t*, and X±{t > t*) = 
a+a-fo*-*^ when t > t*. Essentially, the point of this approximation is that in regions of the 
Brillouin zone where the induced gap is smaller than u, the second sublattice is treated as 
normal, while in the remaining regions, uj is treated as much smaller than the induced gap. 
We find that pi(co') ~ pi(^) + Pi(^) with 



p'lM 



AAqLU 1 



K 



A+-A- 
A4- 



(Al6a) 



and 



P» 



a2tl 



TT' 



'[yUsCTi -/Xi(J2]^ 



+ 



COS — 

2 ^to^ 

4aot*2 



K 


r 21 
1 ^' 




t* 

'To 


M^ 


2 



7r2[/X2Cri - yUl(T2]2 [ 4(72 



I/U2I 



cos-i(-) 



(Al6b) 



when t*{u!) < to, and 



P'l 



(T2tl 



:K 



1- 



j4_ 



(Al6c) 



27r2[/i2cri -/iicr2]^' 

when t*{ij) < to- Equation (|A16a|) is just the DOS which would result in the decoupled 
(to = 0) limit. Equation ( [A16b|) shows that the interlayer coupling induces an inner gap for 
t*(u;) < to (which is the prediction made in Eqn. (pS])). The width of this gap is proportional 
to tg and is difficult to see for weak interlattice coupling. 

We can find the DOS in the second sublattice, p2, in an analogous manner: 



P2(^) 



cr27r-' 



:K 



1- 



P2 



I60-J 



COS — 



'to' 



27 



+ 



At 



*2 



2+2 



a2TT% 



4(T2 






utldi 
47rAo 



[/i20-i - yUiCr2 



-2 



1- 






(A17a) 



when t* < to and 



P2IWJ 



o-27r^ 



ris: 






CUtnO"! r n 

+ -^-^[Ai20-1 -/ilO-2j 



47rA, 







1- 



16af 



(A17b) 



when t* > to- Again we see that there is an induced gap for t*{uj) < to, due to the interactions 
of the sublattices. For t*{uj) > to, P2(^) is just the normal state DOS at the Fermi surface 
with corrections of order tp due to the adjacent layers. In Fig. |^ we graph the low energy 
DOS calculated from equations (p!4D, ( |A10|) and from our analytical approximations. We see 
that, although our analytic expressions for pi and p2 tend to overestimate the size of the 
inner gap, they contain the essential description of the DOS. 
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FIGURES 
FIG. 1. Tc and Ao(T = ^)/Tc vs to- The model parameters are cJi = 1, cr2 = 0.6, /ii = —0.8, 

112 = 0.4 

FIG. 2. Quasiparticle energy dispersion along the line k^ = ky, k^ = 0. We have chosen an 
s-wave gap for purposes of illustration since a d-wave gap vanishes along this line. The features 
which are important for our discussion of the density of states are the two minima of E^ , which 
are related to the intrinsic and induced gaps, and the two avoided band crossings. The parameters 
are ai = 1, o"2 = 0.6, /^i = A*2 = —0.8, to = 0-2, Ao = 0.1 

FIG. 3. Density of states in the normal and superconducting sublattices for the uncoupled limit. 
There are logarithmic divergences at a; ~ — Aii,— /^2 due to the intrinsic van Hove singularities of 
^1 and ^2- There is also a d-wave gap at the Fermi surface in the first sublattice. The choices of 
parameters are ui = 1, cr2 = 0.6, fJ-i = fJ-2 = 0.8, to = 0.01, Ao = 0.1, Tc = 0.086 

FIG. 4. Density of states in the normal and superconducting sublattices. Here the bands ^i 
and ,^2 cross away from their van Hove singularities. In addition to the intrinsic singularities at 
uj ~ —fj,i,—fi2 and the d-wave gap in pi there are two gap-like structures at w ~ ±0.55. There 
is also an induced gap at the Fermi surface in p2 which is too small to show on this plot. The 
parameters are ui = 1, o"2 = 0.6, pi = —0.8, fi2 = 0.4, to = 0.1, Ao = 0.09, Tc = 0.086 

FIG. 5. Density of states in the normal and superconducting sublattices. This figure shows the 
effect of increasing the interlayer coupling on Fig. ^ The smearing of the logarithmic singularities 
by the third dimension is clearly evident. The induced gap at the Fermi surface is now visible 
(although the fact that /?2(0) = is not clear in the plot), and the gap-like structures away from 
the Fermi surface are much larger, and shifted slightly, as we would expect from Eqn. (|2^). Here 
we have ai = 1, 0-2 = 0.6, 1^1 = -0.8, fi2 = 0.4, to = 0.4, Ao = 0.15, T^ = 0.086 
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FIG. 6. Density of states in the normal and superconducting sublattices. In this case, the van 
Hove singularities intrinsic to the sublattice dispersions at cj ~ 0.8 are strongly suppressed by 
an avoided band crossing. This should be compared with Fig. ^, where the interlayer coupling is 
the same, but the intrinsic van Hove singularities of ^i and ^2 can still be identified. The other 
interesting feature of this figure is that we can see clearly that the induced gap is reflected in pi . We 
should also note again that p2 actually vanishes at oj = 0. Here cJi = 1, (T2 = 0.6, /xi = /X2 = —0.8, 
to = 0.4, Ao = 0.15, Tc = 0.086. 

FIG. 7. Low frequency densities of states: (a) Determined numerically from the exact ex- 



pression, Eqn. (14) (b) Determined numerically from the approximate expression, Eqn. ( AlO ) (c 



Determined from the approximate analytical expressions ( |A16| ) and (A17). The model parameters 
are ci = 1, o"2 = 0.6, fJ-i = ^2 = —0.8, to = 0.2, Ao = 0.15. 

FIG. 8. Quasiparticle energy dispersion along the lines (a) kx = ky, k^ = and (b) kx = ky, 
kz = IT for an s-wave gap. In this case, the Fermi surfaces of ^1 and ^2 coincide. Whether or not 
E-, has a single minimum or double minimum structure depends on the value of t(k) (see eqn. 

)). The parameters are ai = 1, (T2 = 0.6, pi = —0.8, ^2 = —0.48, to = 0.2, Ao = 0.1. 



FIG. 9. Density of states. The Fermi surfaces of .^1 and ^2 are coincident. In addition to the 
intrinsic singularities at a; ~ — /xi, -^2; there is a nested gap structure at the Fermi surface. The 
outer gap is related to the usual d-wave gap in a 2-dimensional superconductor while the inner 
gap comes from a mixture of band structure and superconductivity effects. The parameters are 
cTi = 1, CT2 = 0.6, fii = -0.8, ^2 = -0.48, to = 0.1, Ao = 0.09, T^ = 0.086. 

FIG. 10. Density of states. Again the Fermi surfaces of ^i and .^2 are coincident. The increased 
interlayer coupling has changed the appearance of the DOS considerably from Fig. p. The outer 
gap no longer exists and the gap which remains has a complicated dependency on both to and Aq. 
This DOS would not be simple to distinguish from that of a 2-dimensional d-wave superconductor 
in a tunneling experiment, even though the physics of the two systems is different. The parameters 
are ai = 1, 0-2 = 0.6, pi = -0.8, ^2 = -0.48, to = 0.4, Ao = 0.15, T^ = 0.086. 
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FIG. 11. Feynman diagram representation of the equations for Gu: (a) and (b) are diagram- 
matic forms of Eqns. ( |50aD and ( ^Ob| ) respectively. These equations form a self consistent set. In 
(c), the final term in (a) is iterated once, using (b), to give the lowest order effect of the mixing of 
interlayer coupling and the mean-field. Thick lines represent full Green's functions, thin lines rep- 
resent uncoupled Green's functions, rightward pointing arrows represent spin up electrons, leftward 
pointing arrows represent spin down holes and vertical dashed lines represent interlayer hopping. 
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